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Abstract 

The significance of final-state interactions in — > vrvr decays is plienomenologically 
demonstrated by taking elastic tttt ^ tttt and inelastic vrvr ^ DD rescattering effects 
into consideration. We find that the present experimental data on — > vr^vr" can well 
be understood in this approach without fine-tuning of the input parameters, and large 
CP-violating asymmetries are expected to manifest themselves in such charmless rare 
processes. 
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Recently the branching ratio of tt+tt has been measured independently by CLEO, 

BABAR and BELLE Collaborations: 



^4.3l}:^ ± 0.5) X 10-6 (CLE0[1]) , 
B{B^d - ^^^n = { (6.3111 ± 1-6) X 10-6 (BELLE0) , (1) 

[d.Smtli] X 10-6 (BABARi) . 



Theoretical predictions for B{B^ tt+tt-), as those given in Refs. @, ||, |^ based on the 
QCD-improved factorization, are in good agreement with the BABAR data but difficult to 
coincide with the CLEO data. It has to be seen, in the near future, how three measurements 
will reach full consistency and whether the final experimental result can well be understood in 
the factorization approach without fine-tuning of the input parameters. 

If the present CLEO data are taken seriously, it seems necessary to modify the theoretical 
prediction for B{B'^ — tt^tt-). To do so one naturally speculates that final-state interactions in 
S — > TTTT decays might be significant and ought not to be ignored in the factorization approach. 
In Refs. p, 3 the elastic tttt ^ hit rescattering effects on the branching ratios and CP asym- 
metries of i? — > TTTT transitions have been demonstrated to be important. It is also likely that 
such rare nonleptonic processes are deeply involved in inelastic final-state interactions |10|, |Tl|. 
However, it remains a big challenge today to handle the inelastic rescattering effects on B decays 
in a systematic and quantitative way |jl2| . 



In this paper we attempt to follow a purely phenomenological approach to demonstrate the 
infiuence of inelastic final-state interactions on Bd — >■ vrvr decays and CP violation. The essential 
argument is that there may a priori exist a two-step process Bd —>■ DD =^ tttt, arising from 
inelastic vrvr ^ DD rescattering, in addition to the direct decay mode 5^ — ^ vrvr. We find that 
both elastic and inelastic rescattering effects are possible to modify the predictions based on 
the QCD-improved factorization. More precise measurements of 5 — > vrvr decays are expected 
to clarify whether final-state interactions in them are really significant or not. 

It is well known that the final state vr+vr~ or vr^vr" of Bd decays consists of both 1 = and 
1 = 2 isospin configurations, and D^D^ or D^D^ consists of both / = and 1 = 1 isospin 
configurations. Under inelastic rescattering the / = configuration of vrvr can mix with that of 
DD, leading to a two-step decay mode Bd — * DD =^ vrvr. The final states vr^vr° of B^ decays, 
which only have the 1 = 2 isospin configuration, cannot mix with DD. But vr^vr° are possible to 
mix with the final states like p^p°, and vr+vr- (or vr°vr°) could also mix with the final states such 
as K^K~ and K^K^ [^. For simplicity, we assume that the inelastic final-state interactions of 



Bd vrvr decays are dominated by the 1 = channel mixing via vrvr ^ DD scattering, leaving 
the 1 = 2 state of vrvr unmixed with others. In the assumption made above and in the neglect 
of small electroweak penguin contributions, the amplitudes of B^ — vr+vr", B^ tt^tt^, and 
5^ —>■ vr"'"vr° decay modes can be written as |jl3[ [|: 

v4(fi°^vr+vr-) = V2S^^AY + V2{S^''A^o" + So^K^) . 
^ vr°vr°) = 2^2""^r - (^o'^o" + Sq^^o^) , 

^ vr+vr°) = SA^ , (2) 
where denote the "bare" isospin amplitudes of -B tttt, Aq^ stands for the 1 = 



isospin amplitude of Bd DD |jTj], and Sq are the inelastic- rescattering matrix ele- 



^Note that the sign of Aq"^ is here taken to be different from that in Ref. The present choice wiU prove 
convenient when the factorization approximation is apphed to the isospin ampHtudes. 
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ments connecting the unitarized isospin amplitudes to the bare ones |Tl|. Obviously 5*^ 



and = Sq'^ = 1 held, if there were no mixture between the 1 = states of rni and DD. 

The bare isospin amplitudes 2 and Aq^ can be calculated with the help of the effec- 
tive weak Hamiltonian [^] and the QCD-improved factorization [^. After a straightforward 
calculation, we obtain 

= ^ KbVud (2ai - as + 3< + 3<C) + V^^Kd {3al + 3a^C)] T^e''° , 

= ^V:,Vud{ai + a2)T^e''\ 


= ^[V:,V^A< + <^n) + V^V,,ia, + al + alU]TDe''\ (3) 
in which Vub, Vud, Vcb, and Vcd are the quark mixing matrix elements ; cii, (32; 0,4 

and ag''^ are the 

QCD coefficients independent of the renormalization scheme ; and C,d are the factorization 
parameters arising from the transformation of (V — A) ( V + A) currents into (V — A) (V — A) 
ones for the penguin operators and Qq |T^; 6q and 62 are strong (isospin) phases, and 



denote the factorized hadronic matrix elements B ^ nn and B DD decays respectively. 
Under isospin symmetry and C,^ read 

2m! 



{rrib - mu){mu + m^) 



c = ^JI^ . (4) 

{nib - mc){mc + ma) 



In terms of the relevant decay constants and form factors, one gets 



— m. 



Td = ifDF^''{ml)(ml-ml) . (5) 

It is worth remarking that the isospin amplitudes Aq\ and Aq^ have been calculated separately 
in the factorization approximation, hence the contribution of A^^ to A{B(i tttt) need in prin- 
ciple be normalized. Such a treatment is however unnecessary in the approach under discussion, 
because the normalization factor of Aq^ can always be absorbed into the unknown parameter 
Sq^ in Eq. (2). This point will be seen more clearly later on. 

Following Eq. (2) one may write down the similar isospin relations for the amplitudes 
of B^ —>■ 7r+7r~, B^ — > 7r°7r°, and B~ — > 7r~7r° decay modes, whose bare isospin amplitudes 

2 and Aq^ can directly be read off from 2 and A^^ in Eq. (3) with the replacements 
Vub^ud =^ VubV*^ and V*f,Vcd =^ VcbV*^- Of course {A'^'^/A'^'^l = 1 holds in the approximation 
of neglecting the electroweak penguin effect. 

Let us give a brief retrospection of the conventional calculations of 5 —> tttt decays, in which 
final-state interactions are considered only at the quark level |Q, |^, ^ |^. Taking 6q = 62 = 0, 
5*0 = 0, and Sq'^ = 82'^ = 1 (i.e., neglecting both elastic tttt ^ vrvr and inelastic vrvr ^ DD 
rescattering effects), we arrive from Eqs. (2) and (3) at 

A(S° ^ 7r+7r-)o = ^[V:,Vudia^ + al + a^,U + V:,V^ial + alU]T., 



^Note that there is a factor of 1/2 for the phase space of tt^tt^ due to identical final-state particles. 

This factor has already been absorbed into the transition amplitude A{B'^ ttOttO) in Eqs. (2) and (6). 
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A(50^7rV)o = ^[V:,V^Aa2-<-aia-V:,V,,{al + alU]T^, 

^ 7r+7r°)o = ^V:,Vud{ai + a2)T^, (6) 

where the subscript "0" denotes the absence of final-state interactions at the hadron level. Such 
a result has been presented in Refs. 0, |^ . The branching ratios of -B — tttt decays can then be 
computed, under isospin symmetry, by using the formula 



B{B ^ 7r7r)o = , \A{B ^ 7r7r)o|' , (7) 

where denotes the i?-meson lifetime. To obtain the quantitative results for E{B —>■ ttii)^, we 
input the following typical values of the quark mixing parameters: \Vud\ = 0.9735, \Vcd\ = 0.224, 



iKbl = 0.0402, \Vub\/\Va\ = 0.090, and 7 = arg(-V;i,V;dKfeK*J = 65° |T6|. Furthermore 



= 1.6 ps, U = 130.7 MeV, and F^^{ml) ^ F(f^(0) = 0.3 are used [|TB], [T^. Following Beneke 
et al. in Ref. @], we adopt ai = 1.038 + O.OlSz, 03 = 0.082 - O.OSOz, < = -0.029 - 0.015i, 
and 04 = —0.034 — O.OOSi at the scale /i = m^, and neglect the formally power-suppressed QCD 
coefficients Qq''^ in the heavy quark limit. The predictions for the branching ratios of 5 ^ vrvr 
turn out to be f\ 

B{B^a 7r+7r-)o ^ 1.0 x 10"^ , 

BiB''^ ^ 7rV)o ^ 1.4 X 10-^ , 
^0 



B{B+ ^ 7r+7r")o ^ 5.8 x 10"^ ; 



and 



B{B^a n+TT-)o ^ 9.8 X 10"^ , 
^ 7rV)o ^ 1.5 X 10-^ , 
B{B- 7i-7i°)o ^ 5.8 X 10-^ . (9) 

First of all we observe that 7i^n^)o deviates almost a factor of 2 from the CLEO 

data given in Eq. (1), but it is in good agreement with the preliminary BABAR data. At 
present it remains too early to claim any discrepancy between the theoretical prediction and the 
experimental measurements. Secondly, the branching ratios of B^ — > n^n^ and B^ — > vr^vr^ are 
rather sensitive to the values of the QCD coefficients and the weak phase 7. The implication 
of B{B^ 'K^'K^)o/B{B^ 7r°7r°)o ~ 0.1 is quite clear: the direct CP-violating asymmetry 
between these two decay modes might be of 0{1). In addition, B^ — > tt'^tt^ transitions can be 
calculated in a relatively reliable way at the tree level, thus a measurement of their branching 
ratios will test the validity of the QCD-factorization approximation. Actually the preliminary 
CLEO data B{B^ 7r±7r°) = (5.4t|J ± 1.5) x 10-*^ [|], though not yet formally announced, do 
agree very well with the theoretical prediction. 

Now we recalculate the branching ratios of Bd tt+tt^ and Bd tt^tt^ decays by taking the 
final-state rescattering effects into account. We assume that the inelastic rescattering matrix 

§The upper bound on the branching ratio of ir'^n'^ has been reported by the CLEO CoUaboration [Q: 

7r°7r°) < 5.7 x 10~^; and the upper bounds on the branching ratio of Bi^ tt^tt^ have been reported 
by both CLEO and BELLE Collaborations: B{B^ 7r±7r°) < 1.27 x IQ-^ (CLEO 0) and B{B^ 7r±7r°) < 
1.01 X 10-5 (BELLE (I). 
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Figure 1: The branching ratio of tt+tt" predicted in the factorization approximation with 

different input values of the inelastic rescattering parameter k and the isospin phase shift 6. 
The dashed region is favored by the CLEO data M. 



S is approximately diagonal (i.e., 



Ctttt 



1) and its off-diagonal element S^^ is only 



of (9(10 ^) in magnitude. In this reasonable assumption, which should not be far away from 
reality, one can simplify Eq. (2) and arrive at 



TT TT 



Gf r 



[ai + a2)e" 



+ -(2ai-a2) + 



ain + {al + ali^)\]T, 



-(ai + a2)e 



v:,v,d 



a2) - « + <C) 



(10) 



where 5 = 82 — Sq and k = Sq Td/T^,. In obtaining Eq. (10), we have neglected the small 
quantities of (9(|a4'^| ■ and (9(|ag'^| ■ \k\). Indeed we find Td/T^, ^ 3 for fjj = 200 MeV 
and F(f^(0) = 0.7 [0, with the help of Eq. (5). Therefore \k\ ~ C(|^o^|) ~ 0(10-^) 
holds. We observe that A{B^ 7r+7r~)o and A{B^ vr°7r°)o in Eq. (6) can be reproduced 
from Eq. (10) by setting k = and 6 = 60 = 0. The free parameters k and 6 are likely 
to affect the branching ratios of Bd tttt transitions significantly. For illustration we take 
K = — 0.05, —0.05i, 0, +0.05i, +0.05, respectively, to compute B{Bd tttt) as a function of 
S G [0,27r] by using Eq. (7). The numerical results are shown in Figs. 1 and 2. 

Fig. 1 indicates that in the absence of inelastic tttt ^ DD rescattering a good agreement 
between the theoretical value of B^B*^ — ^ vr^vr") and the CLEO data invokes S ~ 0.57r or 
1.57r; i.e., there may exist significant elastic tttt ^ vrvr rescattering. This result is apparently 
consistent with the analyses made in Refs. p, ^ . Taken the final-state vrvr ^ DD rescattering 
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Figure 2: The branching ratio of tt^tt'^ predicted in the factorization approximation with 

different input values of the inelastic rescattering parameter k and the isospin phase shift 5. 



effect into account, the naive value of B{B^ — > 7r+7r~)o given before can be lowered even in the 
case (5 = 0: e.g., we obtain B{B^ — > tt^tt~)/B{B'^ 7r+7r")o ~ 0.7 by taking Ren = and 
ImK = 0.05. We therefore conclude that both kinds of final-state interactions are important 
and non-negligible. Indeed there is a rather large (k, (5)-parameter space, as shown in Fig. 1, in 
which the present data of CLEO, BELLE or BABAR on B^ tt+tt" and B^ tt+tt" decays 
can well be accommodated. 

Fig. 2 shows our prediction for B{B^ —>■ tt^tt^) in the presence of both tttt ^ tttc and 
TTTT ^ DD rescattering effects. We observe that the naive value of -B(-B° —>■ 71^71^)0 can be 
enhanced up to two orders of magnitude. In the region of 5 ~ O.bn or 

2 X 10^^ to 8 X 10"^ is expected for different values of k. Considering the CLEO upper limit 
B{B^ n^Ti^) < 5.7 X 10~^ however, we see that the region 6 < O.Stt or 5 > I.Stt is more 
favored 0. An experimental determination of the branching ratios of B^ tt^tt^ and B^^ vr^vr" 
transitions at the level of 10~^ should signify significant final-state interactions. 

We proceed to calculate the direct CP-violating asymmetries in B^ vs B^ tt+tt^ and 
decays, defined respectively as 



|A(i?g^7r+7r-)p-|A(^°^7r+7r-)p 
|A(i?o^vr+7r-)P + |A(50^7r+vr-)P 



These asymmetries can be observed time-independently on the T{AS) resonance with a trivial 
dilution factor 1/(1 + x^) ~ 0.66 due to B^-B^ mixing [1^, or time-dependently at asymmetric 



^Note that 5 « 11° is expected in the Regge model IH 
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Figure 3: The CP-violating asymmetry of vs 7i'^7t~ decays predicted in the factoriza- 

tion approximation with different input values of the inelastic rescattering parameter k and the 
isospin phase shift 6. 



S-meson factories running around the T(4S') energy threshold 20 . The numerical results of 



.4(7r+7r^) and ^(vr^vr'') are shown in Figs. 3 and 4, where we have used the same values as 
before for the relevant input parameters. Some comments are in order: 

(a) In the absence of final-state interactions at the hadron level (i.e., S = and k = 0), the 
CP asymmetries ^(tt^tt") ^ 3% and A{7t'^7i^) ~ —82% are a consequence of the interference 
between tree-level and penguin amplitudes, where the non-trivial strong phase shift arises from 
the penguin quark-loop function Note that we have neglected possible effects from the elec- 



troweak penguins |2^, the space-like penguins |2^, and the self- interference of different penguin 
loops, as they are generally expected to be insignificant in the transitions under consideration. 

(b) If only the vrvr ^ tttt rescattering effect is "switched on", ^(tt+tt") undergoes an os- 
cillation with increasing values of 6 and its magnitude can be as large as 25% for 6 ~ 1.27r, 
while the magnitude of ^(vr'^Tr^) always decreases when 6 deviates from or 27r. In this case 
CP violation remains resulting from the interference between tree-level and penguin amplitudes, 
but the relevant isospin phase differences may play a more important role than the strong phase 
shifts induced by penguin loops at the quark level. 

(c) If only the vrvr ^ DD rescattering effect is "switched on", the magnitude of ^(vr+vr") 
can remarkably be enhanced (e.g., A{it^7i~) ~ 35% for k = — 0.05i), but that of ^(vr^vr^) 
becomes smaller than in the case k = 0. There are two sources of CP violation: one is the 
interference between tree-level and penguin amplitudes, and the other is the interference between 
two different tree-level amplitudes as a result of tttt ^ DD rescattering. The latter is measured 
by K, whose magnitude and phase can both affect the CP-violating asymmetries in a significant 
way. 

(d) In general nn ^ tttt and tttt ^ DD rescattering effects should both be taken into account. 
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Figure 4: The CP-violating asymmetry of vs n n decays predicted in the factorization 

approximation with different input values of the inelastic rescattering parameter k and the 
isospin phase shift 6. 



It is then possible to have |^(7r+7r^)| ~ 0{1) and |^(7r°7r°)| ~ 0{1) for appropriate values of the 
input parameters. While large CP asymmetries are likely to be measured in B^ — > vttt decays, 
to pin down the true mechanism of CP violation will be very difficult. Furthermore, the indirect 
CP-violating signals in such decay modes (arising from the interplay of direct decay and B^-B^ 



mixing) are unavoidably contaminated by significant final-state rescattering effects |jT3|, |2^. It 
is therefore a big challenge to extract any information on the weak CP-violating phases from 
Bd — >■ TTTT transitions. 

In summary, we have demonstrated the significance of elastic and inelastic final-state rescat- 
tering effects in Bd tttt decays. Our treatment of the complicated inelastic final-state in- 
teractions is so ley to take the simple vrvr ^ DD rescattering into account, hence it remains 
quite preliminary and can only serve for illustration. Nevertheless, the present experimental 
data on 5 — > vrvr decays can well be understood in our approach without fine-tuning of the 
input parameters, and large CP-violating asymmetries are expected to manifest themselves in 
such charmless rare processes. We remark that further effort is desirable towards a deeper 
understanding of the dynamics of nonleptonic B decays. 

I would like to thank X. Calmet, J.N. Ng, X.Y. Pham, and in particular D.S. Yang for useful 
discussions. I am also grateful to J.M. Gerard and W.S. Hou for useful comments. 
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